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We study a dynamically generated pattern in height gra- 
dients, centered around the active growth site, in the steady 
state of a self-organised interface depinning model. The pat- 
tern has a power-law tail and depends on interface slope. An 
approximate integral equation relates the profile to local inter- 
face readjustments and long-ranged jumps of the active site. 
The pattern results in a two-point correlation function satu- 
rating to a finite value which depends on system size. Pattern 
formation is generic to systems in which the dynamics leads 
to correlated motion of the active site. 

PACS numbers: 47.54.+r, 68.10.Gw, 05.40.+j, 47.55.Mh 

A common feature of many self-organised systems is 
that their dynamics is extremal JjJ: that is, activity is 
triggered at that spot where a link is weakest, or a force 
the strongest, and not statistically uniformly over the full 
system. As a result, the active site performs an erratic 
motion with nontrivial correlations in space and time . 
In this Letter, we demonstrate an intriguing and generic 
feature of such systems: the existence of a well-defined 
spatial pattern around this dynamical center of activity. 

We study this phenomenon in a model of a self- 
organizing interface in a random medium. In steady 
state, we find that there is a nonzero time-averaged spa- 
tial pattern of height differences, provided that it is re- 
ferred always to the location of the active site at that 
instant. Figure 1 is a schematic representation of the 
interface profile, time-averaged in a frame which keeps 
the active site at the origin. Interfaces in which r — > — r 
symmetry is present (Fig. la), or absent (Fig. lb), are 
shown, symmetry-breaking being induced by tilting the 
interface ||. In the former case, the profile has a cusp at 
the origin, while in the latter case the dominant feature 
is a larger slope near the origin — a local amplification 
of the broken symmetry. The activity-centered pattern 
(ACP) is defined in terms of height gradients, 

y(r) = ±[(Vh(r + R(t)))-m], (1) 

where R(t) is the position of the active site at time t, 
h(r') denotes the height at the site r', (...) is a time aver- 
age and m is the overall slope of the interface. ^From 
numerical simulations, we find that \&(r) has a long- 
ranged (power law) tail with exponents which differ in 
tilted and untilted interfaces. We also find that the ACP 
leads to an unusual finite size effect, namely the satura- 
tion of a two-point correlation function with increasing 
separation. Finally, we present an approximate integral 
equation which relates the ACP to the local dynamics 



of interface readjustment and the long-ranged dynamics 
of active-site motion. This provides an understanding of 
the form of the induced pattern as a function of spatial 
correlations in the activity. 

In our model the interface is taken to be a directed 
path on a square lattice (Fig. 2), with tilted cylindri- 
cal boundary conditions ||] which ensure that the mean 
slope is preserved. To every bond k on the lattice is 
pre-assigned a fixed random number drawn from the 
interval [0,1]. Motion of the interface is initiated at that 
bond, on the forward perimeter of the interface, which 
carries the smallest random number. In order that the 
length of the interface not change (an effect of surface 
tension), its directed walk character is preserved by local 
re-adjustments; if the chosen minimal bond has a posi- 
tive (negative) slope, the sequence of links with negative 
(positive) slope just below (on the left) also advance ||, 
as illustrated in Fig. 2. This model resembles that in- 
troduced by Sneppen || but differs from it in that the 
length of the interface is a strict constant of the motion 
and the fk's are associated with bonds rather than sites. 

The interface dynamics is equivalent to the dynamics of 
a system of hard-core particles on a ring, with a positive- 
slope link of the interface represented by a particle (rij = 
1) and a negative-slope link by a hole (rij = 0); see Fig. 
2. The difference in height of the interface between sites 
h and j 2 is given by h h - h jl = Vj : ' [2u , - 1). Each 
site j of the ring also carries the random number fj as- 
signed to the bond in front of the corresponding spot on 
the interface. In each time step, activity is initiated at 
the site with the minimum fj . If this site contains a par- 
ticle (hole), it exchanges with the first hole (particle) on 
the left (right). All sites hopped over, including the two 
which exchange the particle and hole, are refreshed by 
assigning a new set of fj's. The overall particle density 
p is strictly conserved, and determines the mean slope 
m = 2p — 1 of the interface. There is a nonzero current, 
which translates into a forward-advancing interface. 

^From time to time, the interface aligns along directed 
spanning paths in a percolation problem, as in the Snep- 
pen model |,|. In our case, the corresponding perco- 
lation problem is simpler, viz. directed bond percola- 
tion (DP) U on the dual lattice, or equivalently, diode- 
resistor percolation (DRP) [[l0| on the original lattice. 
For a given configuration {fk} and a trial threshold value 
fa, occupy all bonds k' with fy > f by diodes (one- 
way connections pointing against the direction of advance 
of the interface, see Fig. 2), and place resistors (two- 
way connections) on the remaining bonds. The diode 
concentration is p = 1 — /„. Infinite connected paths 
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(stoppers) first form along the easy direction of DP at 
P = Pdp — 0.6445 on the square lattice. An untilted 
interface (m = 0) aligns along stoppers in the course of 
its motion For m/0, the corresponding threshold 

p c (m) > pdp is that value of p for which the edge of the 
DP cone § has slope m. As it evolves, a tilted interface 
of slope m aligns along stoppers of that slope . 

The pattern defined in (1) is associated with a spa- 
tial profile of the density in the particle-hole model, as 
the height gradient maps on to the density. The defining 
equation (1) now reads W(r) = (n(r + R(t))) — p. We 
studied \t(r) by Monte Carlo simulation. In the untilted 
case, \I/(r) is an odd function (Fig. 3a) decaying asymp- 
totically as a power law \r\~~ 6 with 9 — 0.90 ± 0.03. In 
the tilted case (Fig. 3b), as there is no r — > — r symme- 
try, \& (r) does not have a definite parity. It is useful to 
separately analyse \l/±(r) = (*(r) ± r))/2. The odd 
part decays as *_(r) — \r\' e - with 0_ = 1.04 ± 0.05. 
The even part ^+(r) « — b(L) + a \r\~ e where 9 + = 
0.46 ± 0.05 and b(L) — > as the lattice size L — ► oo. 

The density profiles of Fig. 3 correspond to the height 
patterns shown in Fig. 1. In the untilted case, on aver- 
age, the active site is located at the peak (Fig. la) where 
/fc's which have not been sampled earlier are most likely 
to occur. In the tilted case there is also a bootstrap effect 
at work. If p > 0.5, the active site is more likely to con- 
tain a particle. Thus regions to the left are more often 
refreshed, making the active site likely to move in this 
direction and find itself amidst a particle cluster. This 
leads to densities higher than p on both sides of the active 
site (Fig. 3b) for p sufficiently different from 0.5. 

More quantitatively, we may write an integral equation 
for \I>(r) for untilted and tilted interfaces. If the pattern 
is centered at R(t) at time i, the dynamics causes two 
changes at the next instant: (i) A short-ranged read- 
justment of the interface changes the profile near R(t). 
As a result, the local density to the right of the active 
site becomes less than the density to its left. The aver- 
age density change at site R(t) + r defines the density- 
increment function <f>(r). (ii) The active site jumps a 
distance I = R(t+ 1) — R(t); there is a substantial proba- 
bility P(l) for large jumps. Since the ACP is centered at 
the active site, the result of (i) followed by (ii) is that the 
average profile reproduces itself, except that it is centered 
at the shifted site R(t) + 1. Both effects are incorporated 
into the integral equation 



P+(l), decays asymptotically as P+(l) 



~ n + with 



*(r) 



+ 00 



[V{r-l) + $(r-l)]P(l)dl. 
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P(l) decays as a power for large I: P(l) ~ \l\~* (Fig. 
4). In the untilted case P(l) — P(—l) holds because 
of r — > — r symmetry. We find it = 2.25 ± 0.05 which 
compares well with earlier determined values of 7r for the 
Sneppen model @^§. In the tilted case (p ^ 1/2), 
P(l) is not a symmetric function (Fig. 4) and it is con- 
venient to separately analyse the even and odd parts 
P± = (P(l) ± P(-l))/2. We find that the even part, 



7r + = 2.00±0.02. The odd part P_ (I) changes sign (as im- 
plied by the crossing of the curves in Fig. 4) and asymp- 
totically follows P-(l) ~ with tt_ = 2.49 ± 0.06. 
We verified that the values of tt + and 7r_ are the same 
for various p ^ 1/2. 

Equation (2) can be solved by Fourier transform. 
Defining = f+°° e 27ri « r *(r)dr etc we find 
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This equation indicates that the ACP arises as a nonlocal 
response to the local density-increment function $(r). 
Notice that the denominator vanishes as q — > 0. We 
have solved (3) numerically, but it is more instructive to 
examine the small-g (large-r) behaviour. 

In the untilted case, the r — > — r symmetry implies 
that P(l) is even, while $(r) and the profile ^(r) are 
both odd. <5(r) is a short-ranged function with finite 
first moment (j>\. Hence, to first order, $>(q) ~ i<piQ- The 
power-law tails in P(l) imply that P(q) ~ 1 — Alqj^ -1 . 
Thus we find \&(r) ~ sgn(r) |r|~' 3_,r ) i.e. the profile is 
an odd function, with a power law tail. 

In the tilted case, the absence of r — > — r symmetry 
implies that none of P(l), $(r) and ^(r) has a definite 
parity. Since $(r) is short-ranged, « ifaq + foq 2 

as q — ► 0. There is no 4>o term, as the elementary step 
of hopping a particle or hole conserves particle number, 
implying J Q(r)dr — 0. The q — * behaviour of P(q) is 
determined by the asymptotic power law decays of the 
even and odd parts P±{l) as \l\ — > oo. Thus we have 
P+(q) ~ 1 — A\q\ 7r +~ 1 . We might have expected P_(<?) ~ 



Bq + C sgn{q) \q\ 



but in fact the mean velocity 



J lP(l)dl of the active site vanishes |l3| implying B = 0. 
Thus the integral equation predicts that to leading order, 



both ^+(r) and ^_(r) decay as powers 



~ 0± , with 



9 + + 2n + — 7r_ = 3 and 0_ + tt + = 3 . The prediction tt + + 
(9_ = 3 compares quite well with the numerical values 
3.15 for the untilted case and 3.04 for ^_(r) in the tilted 
case. For ty + (r), however, the numerically determined 
value of 9+ + 2tt + — 7r_(~ 1.97) deviates more from the 
predicted value 3. The reason behind the discrepancy is 
that (2) implicitly ignores correlations e.g. in successive 
jump lengths of the active site. 

To check this point, we numerically studied a model 
with no correlations in successive jump lengths. At ev- 
ery instant the jump length of the active site was taken 
from a power-law distribution P(l) chosen suitably so as 
to implement the zero velocity constraint. Short-ranged 
adjustments were taken to follow the same particle-hole 
hopping rules as the extremal model. This model is simi- 
lar in spirit to the Levy flight interface model considered 
in p2| . The results in this case were fully consistent with 



the predictions 9 + + 2tt^ 



7r + = 3. Further, 



we also studied cases in which P(l) is not a power law. 
For short ranged P(l), the ACP has the form of a kink 
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function, with a particle (hole) - rich region on the left 
(right) of the active site. In the limit of strictly infinite 
ranged P(l), tantamount to normal stochastic evolution 
H , the ACP ceases to exist. In all cases we compared our 
numerical results with those predicted by (2) and found 
good agreement. 

The ACP has several interesting consequences. For in- 
stance, we expect there to be a larger length of interface 
in a region of fixed size x around the active site, than 
in a region opposite it. Accordingly, we monitored mean 
squared fluctuations of the height around the instanta- 
neous average, in regions around and opposite the active 
site, and found a pronounced difference (factor ~ 2, for 
tilted and untilted cases with x — 256, L = 4096). This 
effect is smaller at stoppers, indicating that the ACP it- 
self is suppressed there. The excess length of interface 
associated with the ACP may provide a useful way to 
identify the active region in experiment. 

Another consequence of the ACP is that being 
a one-point correlation function, albeit an unusual 
one, it has a strong effect on the customary, space- 
time averaged two-point correlation function C(Ar) = 
{(n(r')n(r' + Ar)}} - {(n{r'))(n(r' + Ar))}. Here r' is 
a fixed site on the lattice, (...) stands for a time av- 
erage and {...} stands for an average over all sites r' . 
Numerical results for C(Ar) (Fig. 5) show that it sat- 
urates at a value C sat which decreases with increasing 
size This unusual behaviour can be understood 

in terms of the ACP. Consider the correlation function 
r(r, Ar) = (6n(r + R(t))5n(r + Ar + R(t))) where R(t) 
is the location of the active site , r is the distance 
from the active site and Sn(r + R(t)) = n(r + R(t)) — 
p — ^>(r) is the fluctuation around the average ACP. 
A reasonable expectation is that the <5n's are indepen- 
dent for large separations Ar, i.e. T(r, Ar) — > as 
Ar — » oo. On averaging over r (which has the same 
effect as the average over space-fixed sites r'), we see 



{(p + V(r))(p + #(r + Ar))} 
oo. This predicts the satura- 



that {(n(r)n(r + Ar)}} 
approaches zero as Ar 
tion value C sat to be {(p+^(r))(p+^(r + Ar))} - p 2 . To 
test this, we subtracted this from C(Ar) and found that 
the saturation effect is in fact suppressed strongly (Fig. 
5), supporting our interpretation |I5|. Customarily, the 
large-separation saturation of C(Ar) is associated with 
nonzero values of (n(r')) — p; the unusual aspect here is 
that there is saturation even though (n(r')) = p. 

Though we have focussed on the ACP associated with 
a conserved quantity like the density, a conservation law 
is not essential for the occurrence of a nontrivial ACP. 
For instance, on monitoring the distribution of / r 's, with 
r referred to the active site, we find profiles with power- 
law tails in the interface depinning model with and with- 
out tilt |L6]]. We also find similar patterns in the Bak- 
Sneppen model of evolution p7[ and a modified asym- 
metric version of it jl6| . This indicates that pattern for- 
mation is a generic feature of systems evolving through 
extremal dynamics. 

In summary, we have explored a new aspect of self- 



organised interface depinning, namely the existence of a 
spatial pattern centered around the active site. We have 
shown that the pattern in height gradients has power law 
tails which are sensitive to r — > — r symmetry, broken by 
tilting the interface, and that it explains the saturation of 
the two-point correlation function with increasing separa- 
tion. The integral equation (2) provides an approximate 
description of the formation of the pattern. We have also 
shown that power-law patterns exist in other quantities, 
for instance the distribution of random numbers around 
the minimum, in both the interface and evolution models. 
Our studies suggest that systems which evolve through 
extremal dynamics are likely to exhibit activity-centered 
patterns. 

We thank Gautam Menon and Deepak Dhar for useful 
discussions, and the referee for his comments. 
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Figure Captions 

FIG. 1. Schematic view of averaged interface profile around 
the active site A in (a) untilted (b) tilted interfaces. Tilt 
breaks r — > —r symmetry, an effect that is amplified near A. 

FIG. 2. The bond model. The tilted interface II' advances 
along the extremal perimeter bond A and locally readjusts to 
align along the dashed line. At the next instant, the activity 
moves from A to A'. The corresponding configuration and 
local moves for the particle-hole model are also shown. SS' 
is a stopper, whose perimeter is fully occupied by diodes. 

FIG. 3. Density profiles in the (a) untilted (p = 0.5) and 
(b) tilted (p = 0.75) cases. Inset: The odd part of the profile 
in tilted (open circles) and untilted (filled circles) cases. We 
used L — 16384 and averaged over 10 8 configurations. 

FIG. 4. Monte Carlo results for the probability distribu- 
tion of the jump of the active site for three different densities, 
p = 0.5 (plus sign), p = 0.75 (circles) and p = 0.84375 (tri- 
angles). If p ^ 0.5, there is a left-right asymmetry and the 
asymptotic slope differs from that for p = 0.5. We used L= 
65536 and averaged over 3.10 9 configurations. 

FIG. 5. Density-density correlations for L = 4096 (open 
circles) and for L = 16384 (squares). The saturation is re- 
duced strongly (filled circles) on subtracting the contribution 
of the ACP. We averaged over 10 configurations. 
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